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Abstract. We consider the core algorithmic problems related to verification of systems with respect to three clas¬ 
sical quantitative properties, namely, the mean-payoff property, the ratio property, and the minimum initial credit 
for energy property. The algorithmic problem given a graph and a quantitative property asks to compute the optimal 
value (the infimum value over all traces) from every node of the graph. We consider graphs with constant treewidth, 
and it is well-known that the control-flow graphs of most programs have constant treewidth. Let n denote the number 
of nodes of a graph, m the number of edges (for constant treewidth graphs m = 0{n)) and W the largest absolute 
value of the weights. Our main theoretical results are as follows. First, for constant treewidth graphs we present an 
algorithm that approximates the mean-payoff value within a multiplicative factor of e in time 0{n ■ log(n/e)) and 
linear space, as compared to the classical algorithms that require quadratic time. Second, for the ratio property we 
present an algorithm that for constant treewidth graphs works in time 0{n ■ log(|a • b\)) — 0(n ■ log(n • W)), 
when the output is f, as compared to the previously best known algorithm with running time 0(n^ • log(n • W)). 
Third, for the minimum initial credit problem we show that (i) for general graphs the problem can be solved in 
0{n^ ■ m) time and the associated decision problem can be solved in 0(n • m) time, improving the previous known 
0(n® • m ■ log(n • W)) and 0{n^ ■ m) bounds, respectively; and (ii) for constant treewidth graphs we present 
an algorithm that requires 0(n • logn) time, improving the previous known 0{n* ■ log(n • W)) bound. We have 
implemented some of our algorithms and show that they present a significant speedup on standard benchmarks. 


1 Introduction 


Boolean vs quantitative verification. The traditional view of verification has been qualitative (Boolean) that classifies 
traces of a system as “correct” vs “incorrect”. In the recent years, motivated by applications to analyze resource- 
constrained systems (such as embedded systems), there has been a huge interest to study quantitative properties of 
systems. A quantitative property assigns to each trace of a system a real-number that quantifies how good or bad the 
trace is, instead of classifying it as correct vs incorrect. For example, a Boolean property may require that every request 
is eventually granted, whereas a quantitative property for each trace can measure the average waiting time between 
requests and corresponding grants. 

Variety of results. Given the importance of quantitative verification, the traditional qualitative view of verifica¬ 
tion has been extended in several ways, such as, quantitative languages and quantitative automata for specification 
languages 02811 81 17122116I46I29I : quantitative logics for specification languages imEl; quantitative synthesis for 
robust reactive systems mm ; a framework for quantitative abstraction refinement 03; quantitative analysis of 
infinite-state systems 0231191 : and model measuring (that extends model checking) 03 . to name a few. The core al¬ 
gorithmic question for many of the above studies is a graph algorithmic problem that requires to analyze a graph wrt 
a quantitative property. 

Important quantitative properties. The three quantitative properties that have been studied for their relevance in 
analysis of reactive systems are as follows. First, the mean-payoff property consists of a weight function that assigns to 
every transition an integer-valued weight and assigns to each trace the long-run average of the weights of the transitions 
of the trace. Second, the ratio property consists of two weight functions (one of which is a positive weight function) 
and assigns to each trace the ratio of the two mean-payoff properties (the denominator is wrt the positive function). 
The minimum initial credit for energy property consists of a weight function (like in the mean-payoff property) and 
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assigns to each trace the minimum number to be added such that the partial sum of the weights for every prehx of the 
trace is non-negative. For example, the mean-payoff property is used for average waiting time, worst-case execution 
time analysis l|18|H|19 | ; the ratio property is used in robustness analysis of systems 0; and the minimum initial credit 
for energy for measuring resource consumptions Go). 

Algorithmic problems. Given a graph and a quantitative property, the value of a node is the inhmum value of all traces 
that start at the respective node. The algorithmic problem (namely, the value problem) for analysis of quantitative prop¬ 
erties consists of a graph and a quantitative property, and asks to compute either the exact value or an approximation of 
the value for every node in the graph. The algorithmic problems are at the heart of many applications, such as automata 
emptiness, model measuring, quantitative abstraction rehnement, etc. 

Treewidth of graphs. A very well-known concept in graph theory is the notion of treewidth of a graph, which is a 
measure of how similar a graph is to a tree (a graph has treewidth 1 precisely if it is a tree) ll42l . The treewidth of a 
graph is dehned based on a tree decomposition of the graph ll32l . see Section]^ for a formal dehnition. Beyond the 
mathematical elegance of the treewidth property for graphs, there are many classes of graphs which arise in practice 
and have constant treewidth. The most important example is that the control flow graphs of goto-free programs for 
many programming languages are of constant treewidth ill, and it was also shown in m that typically all Java 
programs have constant treewidth. For many other applications see the surveys lai. The constant treewidth property 
of graphs has also played an important role in logic and verihcation; for example, MSO (Monadic Second Order 
logic) queries can be solved in polynomial time Il25l (also in log-space ||30l) for constant-treewidth graphs; parity 
games on graphs with constant treewidth can be solved in polynomial time and there exist faster algorithms for 
probabilistic models (like Markov decision processes) Ga. Moreover, recently it has been shown that the constant 
treewidth property is also useful for interprocedural analysis GSl- 
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Table 1: Time complexity of existing and our solutions for the minimum mean-cycle value and ratio-cycle value 
problem in constant treewidth weighted graphs with n nodes and largest absolute weight W, when the output is the 
(irreducible) fraction f 0. 
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Table 2: Complexity of the existing and our solution for the minimum initial credit problem on weighted graphs of n 
nodes, m edges, and largest absolute weight W. 


Previous results and our contributions. In this work we consider general graphs and graphs with constant treewidth, 
and the algorithmic problems to compute the exact value or an approximation of the value for every node wrt to 
quantitative properties given as the mean-payoff, the ratio, or the minimum initial credit for energy. We hrst present 
the relevant previous results, and then our contributions. 

Previous results. We consider graphs with n nodes, m edges, and let W denote the largest absolute value of the weights. 
The running time of the algorithms is characterized by the number of arithmetic operations (i.e., each operation takes 
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constant time); and the space is characterized by the maximum number of integers the algorithm stores. The classical 
algorithm for graphs with mean-payoff properties is the minimum mean-cycle problem of Karp ll35l . and the algorithm 
requires 0{n ■ m) running time and 0{n?) space. A different algorithm was proposed in that requires 0(n ■ m) 
running time and 0(n) space. Orlin and Ahuja 1401 gave an algorithm running in time 0{y/n ■ m ■ log(n • W)). For 
some special cases there exist faster approximation algorithms l20l . There is a straightforward reduction of the ratio 
problem to the mean-payoff problem. For computing the exact minimum ratio, the fastest known strongly polynomial 
time algorithm is Burns’ algorithm lT3l running in time 0{n? ■ m). Also, there is an algorithm by Lawler iJTl that 
uses 0(n ■ m ■ log(n • W)) time. Many pseudopolynomial algorithms are known for the problem, with polynomial 
dependency on the numbers appearing in the weight function, see fT7\ . For the minimum initial credit for energy 
problem, the decision problem (i.e., is the energy required for node v at most c?) can be solved in 0{ti? ■ m) time, 
leading to an • m ■ log(n • W)) time algorithm for the minimum initial credit for energy problem ifTOl . All the 
above algorithms are for general graphs (without the constant-treewidth restriction). 


Our contributions. Our main contributions are as follows. 

1. Finding the mean-payoff and ratio values in constant-treewidth graphs. We present two results for constant 
treewidth graphs. First, for the exact computation we present an algorithm that requires 0{n ■ log(|a • &|)) time and 
0{n) space, where | 7 ^ 0 is the (irreducible) ratio/mean-payoff of the output. If | = 0 then the algorithm uses 
0{n) time. Note that log(|a • &|) < 2 log(n • W). We also present a space-efficient version of the algorithm that 
requires only 0(log n) space. Second, we present an algorithm for hnding an e-factor approximation that requires 
0{n ■ log(n/e)) time and 0{n) space, as compared to the ® • log(n • W)) time solution of Orlin & Ahuja, 
and the O(n^) time solution of Karp (see Table[^. 

2. Finding the minimum initial credit in graphs. We present two results. First, we consider the exact computation 
for general graphs, and present (i) an 0 (n • m) time algorithm for the decision problem (improving the previous 
known 0{iff ■ m) bound), and (ii) an 0{iff ■ m) time algorithm to compute value of all nodes (improving the 
previous known 0{n^ ■ m ■ log(n • W)) bound). Finally, we consider the computation of the exact value for graphs 
with constant treewidth and present an algorithm that requires 0{n ■ logn) time (improving the previous known 
0{n^ ■ log(n • W)) bound) (see Table]^. 

3. Experimental results. We have implemented our algorithms for the minimum mean cycle and minimum initial 
credit problems and ran them on standard benchmarks (DaCapo suit IJI for the minimum mean cycle problem, 
and DIMACS challenges m for the minimum initial credit problem). For the minimum mean cycle problem, 
our results show that our algorithm has lower running time than all the classical polynomial-time algorithms. For 
the minimum initial credit problem, our algorithm provides a signihcant speedup over the existing method. Both 
improvements are demonstrated even on graphs of small/medium size. Note that our theoretical improvements 
(better asymptotic bounds) imply improvements for large graphs, and our improvements on medium size graphs 
indicate that our algorithms have practical applicability with small constants. 


Technical contributions. The key technical contributions of our work are as follows: 

1. Mean-payoff and ratio values in constant-treewidth graphs. Given a graph with constant treewidth, let c* be the 
smallest weight of a simple cycle. First, we present a linear-time algorithm that computes c* exactly (if c* > 0) 
or approximate within a polynomial factor (if c* < 0). Then, we show that if the minimum ratio value v* is the 
irreducible fraction |, then v* can be computed by evaluating 0 (log(|a • b\)) inequalities of the form i^* > v. 
Each such inequality is evaluated by computing the smallest weight of a simple cycle in a modified graph. Finally, 
for e-approximating the value v*, we show that 0 (log(n/e)) such inequalities suffice. 

2. Minimum initial credit problem. We show that for general graphs, the decision problem can be solved with two 
applications of Bellman-Ford-type algorithms, and the value problem reduces to finding non-positive cycles in the 
graph, followed by one instance of the single-source shortest path problem. We then show how the invariants of 
the algorithm for the value problem on general graphs can be maintained by a particular graph traversal of the 
tree-decomposition for constant-treewidth graphs. 
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2 Definitions 


Weighted graphs. We consider^n/fe weighted directed graphs G = (V, E, wt, wt') where V is the set of n nodes, 
E CV X y is the edge relation of m edges, wt : iii —)■ Z is a weight function that assigns an integer weight wt(e) to 
each edge e G E, and wt' : E -G N+ is a weight function that assigns strictly positive integer weights. For technical 
simplicity, we assume that there exists at least one outgoing edge from every node. In certain cases where the function 
wt' is irrelevant, we will consider weighted graphs G = {V, E, wt), i.e., without the function wt'. 

Finite and infinite paths. A finite path P = {ui,... ,Uj), is a sequence of nodes Ui G V such that for all 1 < * < j 
we have {ui,Ui+i) G E. The length of P is |P| — j — 1. A single-node path has length 0. The path P is simple 
if there is no node repeated in P, and it is a cycle if j > 1 and ui = uj. The path P is a simple cycle if P is 
a cycle and the sequence (u 2 , ■. - Uj) is a simple path. The functions wt and wt' naturally extend to paths, so that 
the weight of a path P with |P| > 0 wrt the weight functions wt and wt' is wt(P) = wt(ui, and 

wt'(P) = wt'(Mi, Ui+i). The value of P is defined to be wt(P) = For the case where |P| = 0, we 

define wt(P) = 0, and wt(P) is undefined. An infinite path V = {ui,U 2 , ...) of G is an infinite sequence of nodes 
such that every finite prefix P of P is a finite path of G. The functions wt and wt' assign to P a value in Z U {—c», c»}: 
we have wt(P) = wt(ui, Ui+i) and wt'(P) = oo. For a (possibly infinite) path P, we use the notation it G P to 
denote that a node u appears in P, and e G P to denote that an edge e appears in P. Given a set P C V, we denote 
with P 0 B the set of nodes of B that appear in P. Given a finite path Pi and a possibly infinite path P 2 , we denote 
with Pi o P 2 the path resulting from the concatenation of Pi and P 2 . 

Distances and witness paths. For nodes u,v G V, we denote with d{u, v) = infp:„.w„ wt(P) the distance from u to 
V. A finite path P : u v is a witness of the distance d{u, v) if wt(P) = d{u, v). An infinite path P is a witness of 
the distance d{u, v) if the following conditions hold; 

1. d{u, v) = wt(P) = — 00 , and 

2. P starts from it, and v is reachable from every node of P. 

Note that d{u, u) = c» is not witnessed by any path. 

Tree decompositions. A tree-decomposition Tree(G) = (Vt,Et) of G is a tree such that the following conditions 
hold; 

1. Vt = {Bq, ..., Bn'-i '.Vi Bi C V} and Ub gVt ^ (every node is covered). 

2. For all (it, v) G E there exists Bi gVt such that u,v G Bi (every edge is covered). 

3. For all i, j, k such that there is a bag Pj, that appears in the simple path P^ -w Bj in Tree(G), we have P.j n Bj C 
Pfe (every node appears in a contiguous subtree of Tree(G)). 

The sets Bi which are nodes in Vt are called bags. Conventionally, we call Bq the root of Tree(G), and denote with 
Lv(Pi) the level of Bi in Tree(G), with Lv(Po) = 0. We say that Tree(G) is balanced if the maximum level is 
maxp. Lv(Pi) = 0(log n'), and it is binary if every bag has at most two children bags. A bag P is called the root bag 
of a node ii if P is the smallest-level bag that contains u, and we often use P^ to refer to the root bag of u. The width 
of a tree-decomposition Tree(G) is the size of the largest bag minus 1. The treewidth of G is the smallest width among 
the widths of all possible tree decompositions of G. The following lemma gives a fundamental structural property of 
tree-decompositions. 

Lemma 1. Consider a graph G = {V, E), a binary tree-decomposition T = Tree(G) and a bag B ofT. Denote with 
{Ci)i<i<s the components ofT created by removing B from T, and let Vi be the set of nodes that appear in bags of 

component Ci. For every i f j, nodes u G Vi, v G Vj and P : u v, we have that P C B f fb (i.e., all paths between 

u and V go through some node in B). 

Theorem 1. For every graph G with n nodes and constant treewidth, a balanced binary tree-decomposition Tree(G) 
of constant width and 0(n) bags can be constructed in (1) 0{n) time and space (2) deterministic logspace (and 
hence polynomial time) m- 
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In the sequel we consider only balanced and binary tree-decompositions of constant width and n' = 0{n) bags (and 
hence of height 0(log n)). Additionally, we consider that every bag is the root bag of at most one node. Obtaining this 
last property is straightforward, simply by replacing each bag B which is the root of A: > 1 nodes xi,.. .Xk with a 
chain of bags Bi,... ,Bk = B, where each Bi is the parent of and = BiU {a;i_|_i}. Note that this keeps 
the tree binary and increases its height by at most a constant factor, hence the resulting tree is also balanced. 

Throughout the paper, we follow the convention that the maximum and minimum of the empty set is —oo and oo 
respectively, i.e., max(0) = —oo and min(0) = oo. Time complexity is measured in number of arithmetic and 
logical operations, and space complexity is measured in number of machine words. Given a graph G, we denote with 
T(G) and S{G) the time and space required for constructing a balanced, binary tree-decomposition Tree(G). We are 
interested in the following problems. 

The minimum mean cycle problem ll3^ . Given a weighted directed graph G = (V, E, wt), the minimum mean cycle 
problem asks to determine for each node u the mean value n*{u) = mincGC„ , where is the set of simple 
cycles reachable from u in G. A cycle G with is called a minimum mean cycle of u. For 0 < e < 1, 

we say that a value /i is an e-approximation of the mean value n*{u) if |/r — < e • |^*(m)|. 

The minimum ratio cycle problem 1(331 . Given a weighted directed graph G = (y,E,wt, wt'), the minimum ratio 
cycle problem asks to determine for each node u the ratio value v*{u) = mincgc„ wt(G), where wt(G) = 
and C„ is the set of simple cycles reachable from u in G. A cycle G with wt(G) = ly* is called a minimum ratio cycle 
of u. The minimum mean cycle problem follows as a special case of the minimum ratio cycle problem for wt'(e) = 1 
for each edge e G E. 

The minimum initial credit problem lIlOl . Given a weighted directed graph G = (V, E,wt), the minimum initial 
credit value problem asks to determine for each node u the smallest energy value E(u) G N U {c»} with the following 
property: there exists an inhnite path V = {ui,U 2 ■ ■ ■) with u = Ui, such that for every hnite prehx P of 7^ we have 
E(m) + wt(P) > 0. Conventionally, we let E(m) = c» if no hnite value exists. The associated decision problem asks 
given a node u and an initial credit c G N whether E(m) < c. 


3 Minimum Cycle 


In the current section we deal with a related graph problem, namely the detection of a minimum-weight simple cycle 
of a graph. In Section we use solutions to the minimum cycle problem to obtain the minimum ratio and minimum 
mean values of a graph. 

The minimum cycle problem. Given a weighted graph G = (C, E, wt), the minimum cycle problem asks to deter¬ 
mine the weight c* of a minimum-weight simple cycle in G, i.e., c* = min^gc wt(G), where C is the set of simple 
cycles in G. 

We describe the algorithm IVlinCycle that operates on a tree-decomposition Tree(G) of an input graph G, and has the 
following properties. 

1. If G has no negative cycles, then IVlinCycle returns the weight c* of a minimum-weight cycle in G. 

2. If G has negative cycles, then MinCycle returns a value that is at most a polynomial (in n) factor smaller than c*. 

U-shaped paths. Following the recent work of lfT9l . we dehne the important notion of U-shaped paths in a tree- 
decomposition Tree(G). Given a bag B and nodes u,v G B, we say that a path P : u v is U-shaped in B, if one 
of the following conditions hold: 

1. Either |P| > 1 and for all intermediate nodes w G P, we have B is an ancestor of B^, 

2. or |P| < 1 and B is B^ or By (i.e., B is the root bag of either u or v). 

Informally, given a bag B, a U-shaped path in P is a path that traverses intermediate nodes that exist only in the 
subtree of Tree(G) rooted in B. The following remark follows from the dehnition of tree-decompositions, and states 
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that every simple cycle C can be seen as a U-shaped path P from the smallest-level node of C to itself. Consequently, 
we can determine the value c* by only considering U-shaped paths in Tree(G). 

Remark 1. Let C = {ui,...,Uk) be a simple cycle in G, and uj = argmin„.gc'Lv(ui). Then P = 

{uj,Uj+i,... Uk, ui, ■ • ■, Uj) is a U-shaped path in and wt(P) = wt(G). 

Informal description of MinCycle. Based on U-shaped paths, the work of llT9l presented a method for computing 
algebraic path properties on tree-decompositions with constant width, where the weights of the edges come from a 
general semiring. Note that integer-valued weights are a special case of the tropical semiring. Our algorithm MinCycle 
is similar to the algorithm Preprocess from lfT9l . It consists of a depth-first traversal of Tree(G), and for each examined 
bag B computes a local distance map LDg : B x B ^ ZU{oo} such that for each u,v G B, we have (i) LDb(u., v) = 
wt{P) for some path P : u v, and (ii) LDg < minp wt(P), where P are taken to be simple u v paths (or 
simple cycles) that are U-shaped in B. This is achieved by traversing Tree(G) in post-order, and for each root bag B^ 
of a node x, we update LDp^ (u, v) with LDp^ (u, x) -f LDp^ (x, v) (i.e., we do path-shortening from node u to node 
V, by considering paths that go through x). See Figure [^for an illustration. 

In the end, MinCycle returns minj^ LDp^ {x, x), i.e., the weight of the smallest-weight U-shaped (not necessarily sim¬ 
ple) cycle it has discovered. Algorithm [T] gives MinCycle in pseudocode. For brevity, in line we consider that if 
{m, v} ^ E or {u, u} % Bi for some child Bi of B, then LDp^ (u, v) = oo. 


Algorithm 1: MinCycle 

Input: A weighted graph G = {V, E, wt) and a balanced binary tree-decomposition Tree(G) 
Output: A value c 

1 Assign c oo 

2 Apply a post-order traversal on Tree(G), and examine each bag B with children Bi , B 2 

3 begin 

4 I foreach u,v £ B do 

s I Assign LDfl(u,u) t—min(LDBj(M,i!), LDB 2 (u,u),wt(u,u)) 

6 end 

7 Discard LDsj^, LDbj 

8 if B is the root bag of a node x then 

9 I foreach u, w G B do 

10 I Assign LD'b(u, u) t—min(LDB(u, u), LDb(u, x)-I-LDB(a;, w)) 

11 end 

12 Assign LDfl t— LD^ 

13 Assign c t— min(c, LDbC*, a;)) 

14 

15 end 

16 return c 


In essence, MinCycle performs repeated summarizations of paths in G. The following lemma follows easily from lfT9l 
Lemma 2], and states that LDp (u, v) is upper bounded by the smallest weight of a U-shaped simple u v path in B. 

Lemma 2 ( II191 Lemma 2]). For every examined bag B and nodes u,v G B, we have 

1 . LDb(?i, v) = wt{P) for some path P : u v (and LDb(m, v) = 00 if no such P exists), 

2. LDb(m, v) < minp:u^ii wt(P) where P ranges over U-shaped simple paths and simple cycles in B. 


At the end of the computation, the returned value c is the weight of a (generally non-simple) cycle G, captured as a 


U-shaped path on its smallest-level node. The cycle G can be recovered by tracing backwards the updates of line 10 
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Fig. 1; Path shortening in line 10 of MinCycle. When is examined, LDb^{u,v) is updated with the weight of the 
U-shaped path P = Pi o P 2 . The paths Pi and P 2 are U-shaped paths in the children bags Bi and B 2 , and we have 

LDBiiu,x) = wt(Pi). 


performed by the algorithm, starting from the node x that performed the last update in line 13 Hence, if C traverses k 
distinct edges, we can write 


k 

c = wt(C') = ^ ki ■ wt(ei) 


( 1 ) 


where each is a distinct edge, and ki is the number of times it appears in C. 
Lemma 3. Let h be the height o/Tree(G'). For every ki in Eq. Q, we have ki < 2^. 


Proof. Note that the edge = {ui,Vi) is first considered by MinCycle in the root bag Bi of node Xi, where Xi = 
argmaxj,^g{„. Lv(j/i) (line 10 1 . As MinCycle backtracks from Bi to the root of Tree(G'), the edge Ci can be 
traversed at most twice as many times in each step (because of line 10 once for each term of the sum LDb(m, x) + 
10 ^( 2 ;, v)). Hence, this doubling will occur at most h times, and ki < 2^. □ 


Lemma 4. Let c be the value returned by MinCycle, h be the height o/Tree(G), and c* = mincwt(G) over all 
simple cycles C in G. The following assertions hold: 

1. If G has no negative cycles, then c = c*. 

2. If G has a negative cycle, then 

(a) c < c*. 

(b) |c| = 0 (|c*| ■n-2’^). 


Proof By Remark[^ we have that c* = wt(P) for a U-shaped path P : x x. By Lemma|^ after MinCycle examines 
Bx, it will be c < LDb„, {x, x) < c*, with the equalities holding if there are no negative cycles in G (by the definition 
of c*, as then LDs^ (a;, x) is witnessed by a simple cycle). By line 10 c can only decrease afterwards, and again by the 
dehnition of c* this can only happen if there are negative cycles in G. This proves items and 2a and the remaining 
of the proof focuses on showing that |c| = O (|c*| • n ■ 2^). 


By rearranging the sum of Eq. Q, we can decompose the obtained cycle G into a set of non-negative cycles Gf, 
and a set of k'~ negative cycles G~, and each cycle Gi^ and G~ appears with multiplicity kf and k~ respectively. 
Then we have 


c| = |wt(G)| = 

k'+ k'- 

kt ■ wt(G+ ) 

< 

k'- 

k~ 

k''' k 




< ^ fc- . |wt(G-)| < |c*| • ^ k- <\c*\-Y,k-=0 (|c*| • n . 2 ^) ( 2 ) 
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The first inequality follows from c < 0, the third inequality holds by the definition of c* , and the last inequality holds 
since k'~ < k. Finally, we have ^7 — O {n ■ 2^), since k = 0{n), and by Lemma|^we have k~ < 2^. □ 

Next we discuss the time and space complexity of MinCycle. 

Lemma 5. Let h be the height o/Tree(G'). MinCycle accesses each bag o/Tree(G) a constant number of times, and 
uses 0{h) additional space. 

Proof. MinCycle accesses each bag a constant number of times, as it performs a post-order traversal on Tree(G) 
(line 1^. Because it computes the local distances in a postorder manner, the number of local distance maps LD^ 
it remembers is bounded by the height h of Tree(G). Since Tree(G) has constant width, LD^ requires a constant 
number of words for storing a constant number of nodes and weights in each B. Hence the total space usage is 0{h), 
and the result follows. □ 

The following theorem summarizes the results of this section. 

Theorem 2. Let G = {V,E,wt) be a weighted graph of n nodes with constant treewidth, and a balanced, binary 
tree-decomposition Tree(G) ofG be given. Let c*, be the smallest weight of a simple cycle in G. Algorithm MinCycle 
uses 0(n) time and O(logn) additional space, and returns a value c such that: 

1. If G has no negative cycles, then c = c*. 

2. If G has a negative cycle, then 

(a) c < c*. 

(b) |c| = |c*| • 

4 The Minimum Ratio and Mean Cycle Problem s 


In the current section we present algorithms for solving the minimum ratio and mean cycle problems for weighted 
graphs G = (V, E, wt, wt') of constant treewidth. 

Remark!. If G is not strongly connected we can compute its strongly connected components in linear time 1431 . 
and use the algorithms of this section to compute the minimum cycle ratio v* in every component Qi separately. 
Afterwards, we compute {u) for every node u by iteratively (i) finding the nodes u that can reach the component Qj 
where j = argminj^-*, (ii) assigning i'*{u) = u*, and (iii) removing Qj and repeating. Since these operations require 
linear time, they do not impact the time complexity. 

In light of Remark]^ we consider graphs that are strongly connected, and hence it follows that is the same for 

every node u, and thus we will speak about the minimum ratio ly* and mean fj,* values of G. 

Claim 1. Let u* be the ratio value of G. Then v* > v iff for every cycle G of G we have wt^(G) > 0, where 
wty (e) = wt(e) — wt'(e) • v for each edge e G E. 


Proof. Indeed, for any cycle G we have wt(G) > v* > v. Then 


wt(G) > V 


wt(G) — 1 / > 0 


wt(G) — i' ■ wt'(G) 


> 0 


wt'(G) 

vt(G) — 1 / ■ wt'(G) > 0 ^(wt(e) — wt'(e) ■ u) > 0 ■^=>- wtj,(G) > 0 

eeC 


with the equality holding iff wt(G) = u. 


□ 
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Hence, given a tree-decomposition Tree(G), for any guess v of the ratio value v*, we can evaluate whether v* > v 
by constructing the weight function wt^ —\Nt — v and executing algorithm MinCycle on input Gi, = (V, E, wti/). By 
Item 2a of Theorem|^and Claim[^we have that the returned value c of MinCycle is c > 0 iff wt^{C) > 0 for all 
cycles C, iff v* >v (and in fact c = 0 iff = v). 


Lemma 6. Let G = {V, G, wt, wt') be a weighted graph of n nodes with constant treewidth and minimum ratio 
value V*. Let Tree(G) be a given balanced, binary tree-decomposition of G of constant width. For any rational v, the 
decision problem of whether v* > v (or v* = v) can be solved in 0{n) time and O(logn) extra space. 

Proof. By Claim we can test whether v* > v hy testing whether Gjy = (y,G, wtj^) has a negative cycle. By 
Theorem|^ a negative cycle in Gi, can be detected in 0(n) time and using O(logn) space. □ 


4.1 Exact solution 

We now describe the method for determining the value i/* of G exactly. This is done by making various guesses v 
such that u* > 1 / and testing for negative cycles in the graph Gi, = (V, E, vjt,y). We hrst determine whether v* = 0, 
using Lemma|^ In the remaining of this section we assume that ly* f- 0. 

Solution overview. Consider that v* > 0. First, we either hnd that v* G (0,1) (hence [i/*\ = 0), or perform an 
exponential search of 0(log;/*) iterations to determine j G N~^ such that z/* G [2-’“^, 2-^]. In the latter case, we 
perform a binary search of 0(logiy*) iterations in the interval [2-^“^, 2-^] to determine [z^*J (see Figure]^. Then, we 
can write u* = [u*\ + x, where a: < 1 is an irreducible fraction y. It has been shown im that such X can be 
determined by evaluating 0(log b) inequalities of the form x>v. The case for < 0 is handled similarly. 

Lemma 7. Let v* 0 be the ratio value ofG. The value [z^*J can be obtained by evaluating 0(log \v*\) inequalities 
of the form u* > v. 

Proof. First determine whether v* > 0, and assume w.l.o.g. that this is the case (the process is similar if v * < 0). 
Perform an exponential search on the interval (0, 2 • [z/*J) by a sequence of evaluations of the inequality v* >Vi = 2*. 
After log[z^*J + 1 steps we either have [i/*\ G (0,1), or have determined a j > 0 such that z/* G [z/j_i,z/j]. 
Then, perform a binary search in the interval [z^^-i, uf, until the running interval [£, r] has length at most 1. Since 
i/j — Vj-i = Vj-i < V*, this will happen after at most logfz^*] steps. Then either [z/*J = [£J or [z^*J = [rj, which 
can be determined by evaluating the inequality v* > [rJ. A similar process can be carried out when V* < 0. Figure]^ 
shows an illustration of the search. □ 



vq Vi 


V2 ■■■ Vj-l y* 


^0 



Fig. 2; Exponential search followed by a binary search to determine [v*\ 


Let Tniax = maxe wt'(e) be the largest weight of an edge wrt wt'. Since v* is a number with denominator at most 
(n — 1) • Tmax, it can be determined exactly by carrying the binary search of Lemma |^until the length of the running 
interval becomes at most — p (thus containing a unique rational with denominator at most (n — 1) • Tmax)- 

Then v* can be obtained by using continued fractions, e.g. as in ll3^ . We rely in the work of Papadimitiiou BTlI to 
obtain a tighter bound. 
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Lemma 8. Let v* ^ Q be the ratio value of G, such that v* is the irreducible fraction | € (~1) !)• Then v* can be 
determined by evaluating 0{\ogb) inequalities of the form v* > v. 

Proof Consider that ly* > 0 (the proof is similar when i/* < 0). It is shown in Bdl that a rational with denominator 
at most b can be determined by evaluating 0(log b) inequalities of the form i/* > v. We remark that b is not required 
to be known, although the work of ED assumes that a bound on the denominator of i/* is known in advance. □ 

Theorem 3. Let G = (V, E, wt, wt') be a weighted graph of n nodes with constant treewidth, and X = max„ |a„ • 
such that V* (u) is the irreducible fraction Let f (G) and S (G) denote the required time and space for constructing 

a balanced binary tree-decomposition Tree(G) o/G with constant width. The minimum ratio cycle problem for G can 
be computed in 

1. 0{T{G) + n ■ log(A)) time and 0{S{G) + n) space; and 

2. o\s{G) + logn) space. 

Proof. In view of Remark 2 the graph G is strongly connected and has a minimum ratio value v*. Let v* = \ v*\ + \ 

/ ' M 

with Itl < 1 - By Lemma 7 Ml/* J can be determined by evaluating 0(log|i/*|) = 0(log|a|) inequalities of the 
form V* > V, and by Lemma 8 ^ can be determined by evaluating 0{b) such inequalities. A balanced binary tree- 
decomposition Tree(G) can be constructed once in T(G) time and 5(G) space, and stored in 0{n) space. Tree(G) 
is also a tree-decomposition of every Gjy required by Claim[D By Theorem|^a negative cycle in Gi, can be detected 
in 0(n) time and using O(logn) space. This concludes Item[^ Item is obtained by the same process, but with 
re-computing Tree(G) every time MinCycle traverses from a bag to a neighbor (thus not storing Tree(G) explicitly). 

□ 


Using Theorem[2we obtain from Theorem[^the following corollary. 

Corollary 1. Let G = {V, E, wt, wt') be a weighted graph ofn nodes with constant treewidth, and A = max„ |a„ • 6 „| 
such that v*{u) is the irreducible fraction The minimum ratio value problem for G can be computed in 

1. 0{n ■ log(A)) time and 0(n) space; and 

2. O(logn) space. 

By setting wt'(e) = 1 for each e S £' in Corollary [^we obtain the following corollary for the minimum mean cycle. 

Corollary 2. Let G = (V,E,wt) be a weighted graph ofn nodes with constant treewidth, and A = max„ \p,*{u)\. 
The minimum mean value problem for G can be computed in 

1. 0{n ■ log(A)) time and 0(n) space; and 

2. Oilogn) space. 


4.2 Approximating the minimum mean cycle 

We now focus on the minimum mean cycle problem, and present algorithms for e-approximating the mean value p* 
of G for any 0 < e < 1 in 0{n ■ log(n/e)) time, i.e., independent of p*. 

Approximate solution in the absence of negative cycles. We first consider graphs G that do not have negative cycles. 
Let G be a minimum mean value cycle, and G' a minimum weight simple cycle in G, and note that p* G [0, wt(G')]. 
Additionally, we have 


wt(G') < wt(G) 


wt(C') < ^ 


vt(G) wt(G') < (n) 
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Consider a binary search in the interval [0, wt(C")], which in step i approximates fi* by the right endpoint ni of its 
current interval. The etTor is bounded by the length of the interval, hence fii — fi* < wt(C") • < {n — 1) ■ fi* ■ 2“*. 

To approximate within a factor e we require 

2“* • (n - 1) < e i > log(n) + log(l/e) (3) 


steps. 

Remark 3. Note that for the minimum ratio value we have wt(C") < W ■ n ■ v*, where W = maxeg e wt'(e). For 
e-approximating v* we would need i > log(n • W jt) steps. 


Approximate solution in the presence of negative cycles. We now turn our attention to e-approximating p* in the 
presence of negative cycles in G. Note that uniformly increasing the weight of each edge so that no negative edges 
exist does not suffice, as the error can be of order e • \ W~ \ rather than e • /i*, where W~ is the minimum edge weight. 


Instead, let c be the value returned by Min Cycle on input G. Item 2a of Theorem guarantees that for the weight 
function wt_|c|(e) = wt(e) + |c|, the graph C?-|c| = wt_|c|) has no negative cycles (although it might still 

have negative edges). The following lemma states that fj,* can be e-approximated by e'-approximating the value /i'* of 
G_|c|, for some e' polynomially (in n) smaller than e. 


Lemma 9. Let fi* and fi'* be the value of G and G_\c\ respectively, and e some desired approximation factor of p,*, 
with 0 < e < 1. There exists an f such that if p' is an e'-approximation of p'* in G_|c|, then p = p' — |c| 

is an e-approximation of p* in G. 


Proof. By construction, we have p'* = p* -\- IcI, where c defined above is the value returned by MinCycle on G. Let 

: |c*|-710(1). Note 


2 b 


we have that 


be the weight of a minimum-weight simple cycle in G. By Theoremj^Item 
that |c*| < (n — 1) • |p.*|, hence p'* = p* -\-\c*\ ■ nO(i) < \p*\- a for a = . Let e! = e/a. By e'-approximating 

p!* by p! we have 


y - < e' • \P\ \{d’ - Ic|) - {p'* - |c|)| < e' • \p'*\ ^ \p-p*\<e'-\p*\.a<e-\p*\ 

The desired result follows. □ 

Theorem 4. Let G = (V, E, wt) be a weighted graph of n nodes with constant treewidth. For any 0 < e < 1, the 
minimum mean value problem can be e-approximated in 0{n ■ log(n/e)) time and 0{n) space. 

Proof. In view of Remark[^the graph G is strongly connected and has a minimum mean value p*. First, we construct 
a balanced binary tree-decomposition Tree(G) of G in 0{n ■ logn) time and 0{n) space TheoremLet c be the 
value returned by MinCycle on the input graph G. If c > 0, by Lemma|^we have p* > 0, and by Eq. 0 p* can be 
e-approximated in 0(log(n/e)) steps. If c < 0, we construct the graph G_\c\ = E, wt_|j,|). By Lemmal^ p* can 
be e-approximated by e! approximating the mean value p!* of G_|c|, where e' = Ym- construction, (dZ\c\ does 
not contain negative cycles, thus p'* > 0, and by Eq. (j^ p'* can be approximated in 0(log(n/e')) = 0(log(n/e)) 
steps. By Lemma|^ each step requires 0{n) time. The statement follows. □ 


5 The Minimum Initial Credit Problem 


In the curi'ent section we present algorithms for solving the minimum initial credit problem on weighted graphs G = 
(V, E, wt). We first deal with arbitrary graphs, and provide (i) an 0{n ■ m) algorithm for the decision problem, and 
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(ii) an 0{n^ ■ m) for the value problem, improving the previously best upper bounds. Afterwards we adapt our 
approach on graphs of constant treewidth to obtain an 0 (n • log n) algorithm for the value problem. 

Non-positive minimum initial credit. For technical convenience we focus on a variant of the minimum initial credit 
problem, where energies are non-positive, and the goal is to keep partial sums of path prefixes non-positive. Formally, 
given a weighted graph G = {V, E, wt), the non-positive minimum initial credit value problem asks to determine 
for each node u G V the largest energy value E(t 6 ) € Z<o U {— 00 } with the following property: there exists an 
infinite path V = {ui,U 2 ...) with u = ui, such that for every finite prefix P of V we have E(m) -f wt(P) < 0. 
Conventionally, we let E(m) = —00 if no finite such value exists. The associated decision problem asks given a node 
u and an initial credit c € Z<o whether E(w) > c. Hence, here minimality is wrt the absolute value of the energy. A 
solution to the standard minimum initial credit problem can be obtained by inverting the sign of each edge weight and 
solving the non-positive minimum initial credit problem in the resulting graph. 

We start with some definitions and claims that will give the intuition for the algorithms to follow. First, we define the 
minimum initial credit of a pair of nodes u, v, which is the energy to reach v from u (i.e., the energy is wrt a finite 
path). 

Finite minimum initial credit. For nodes u,v G V, we denote with E^(m) g Z<o U {— 00 } the largest value with the 
following property: there exists a path P : u v such that for every prefix P' of P we have E„(u) + wt(P') < 0. 
Note that for every pair of nodes u,v G V, we have E(t 6 ) > E„(t 6 ) 4- E(r;). Conventionally, we let E„(t 6 ) = —00 if no 
such value exists (i.e., there is no path u v). 

Remark 4. For any u G V,let P : u vhea witness path for Eviu) > — 00 . Then 

E^(u) + wt(P) < 0 Et,(M) < —wt(P) < —d{u,v) 

i.e., the energy to reach v from u is upper bounded by minus the distance from u to v. 

Highest-energy nodes. Given a (possibly infinite) path P with wt(P) < 00 , we say that a node a; G P is a highest- 
energy node of P if there exists a highest-energy prefix Pi of P ending in x such that for any prefix P 2 of P we 
have wt(Pi) > wt(P 2 ). Note that since the weights are integers, for every pair of paths P{, P 2 , it is either |wt(P() — 
wt(P 2 )| = 0 or |wt(P{) — wt(P 2 )| > 1. Therefore the set {wt{Pi)}i of weights of prefixes of P has a maximum, and 
thus a highest-energy node always exists when wt(P) < 00 . The following properties are easy to verify: 

1. If a; is a highest-energy node in a path P : u v, then E„(a;) = 0. 

2. If a; is a highest-energy node in an infinite path P, then E(a;) = 0. 

The following claim states that the energy E(u) of a node u is the maximum energy E^(u) to reach a 0-energy node v. 


Claim 2. For every u GV,we have E(u) = maXt,.E(„)=o ^v{u). 

Proof. The direction E(m) > max„.E(„)^g E„(t 6 ) is straightforward. For the other direction, consider that E(m) > —00 
(trivially, —00 < max^.E(„)^Q E^(m)) and let P be a witness path for E(m). Since E(t 6 ) > — 00 , we have wt(P) < 00 , 
and P has some highest-energy node x, thus E(a;) = 0. Since x is on the witness P of E(w), we have E(m) < Ea;(t 6 ) < 
max„.E(^)=o Et,(M). The result follows. □ 


5.1 The decision problem for general graphs 

Here we address the decision problem, namely, given some node u G V and an initial credit c G Z<o, determine 
whether E(u) > c. The following claim states that if E(m) > c, then a non-positive cycle can be reached from u with 
initial credit c, by paths of length less than n. 
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Claim 3. For every u G V and c S Z<o, we have that E(u) > c iff there exists a simple cycle C such that (i) wt(C') < 
0 and (ii) for every v G C we have that £„(«) > c, which is witnessed by a path Py : u v with |P„| < n. 

Proof. For the one direction, if wt(C') < 0 we have wt(C'‘^) < oo, thus C contains a 0-energy node w. By Claim|^ 
E(m) = E„(it) > Eu,(u) > c. For the other direction, let 7^ be a witness path for E(u), and we can assume 

w.l.o.g. that V does not contain positive cycles. Then for every prefix Py : u v of V we have E(u) + wt{Py) < 0, 
thus E„(rt) > E(u) > c, and the n-th such prefix contains a non-positive cycle C. The result follows. □ 


Algorithm DecisionEnergy. Claim [^suggests a way to decide whether E(m) > c. First, we start with energy c from 
u, and perform a sequence of n — 1 relaxation steps, similar to the Bellman-Ford algorithm, to discover the set Vf 
of nodes that can be reached from u with initial credit c by a path of length at most n — 1. Afterwards, we perform 
a Bellman-Ford computation on the subgraph G \ Vf induced by the set Vf. By Claim we have that E(u) > c 
iff G \ Vf contains a non-positive cycle. Algorithm|2|(DecisionEnergy) gives a formal description. The/or loop in 
lines 6p3 is similar to the procedure ROUND from the algorithm of lITOl . 


Detecting non-positive cycles. It is known that the Bellman-Ford algorithm can detect negative cycles. To detect non¬ 
positive cycles in a graph G with n nodes and weight function wt, we execute Bellman-Ford on G with a slightly 
modified weight function wt' for which wt'(e) = wt(e) — -. Then for any simple cycle C in G we have wt(G) < 0 
iff wt' (G) < 0. Indeed, 


wt'(G) < 0 


E - E - < 0 


eec 


eec 


».(C) < H 

n 


wt(G) < 0 


since |G| < n and wt(G) G Z. 


Algorithm 2: DecisionEnergy 

Input: A weighted graph G = {V, E, wt), a node u G U, an initial energy c G Z<o 
Output: True iff E(w) > c 


// Initialization 

1 foreach u e U do 

2 I Assign D{s) oo 

3 end 

4 Assign D{u) 4— c 

5 Assign Vf t— {u} 

// n —1 relaxation steps to discover Vf 

6 for i t— 1 to n — 1 do 

foreach {v, w) G E do 

if D{w) > D{v) + wt(u, w) and D{v) + wi{v, w) < 0 then 

9 I Assign D(w) D{y) -f wt(u, w) 

10 Assign Uf t— U/ U {w} 

11 

12 end 

13 end 

14 Execute Bellman-Ford onG \ Vf 

15 return True iff a non-positive cycle is discovered 


The correctness of DecisionEnergy follows directly from Claimj^ The time complexity is 0{n ■ m) time spent in the 
for loop of lines 6p3 plus 0{n ■ m) time for the Bellman-Ford. We thus obtain the following theorem. 
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Theorem 5. Let G = (V, E, wt) be a weighted graph of n nodes and m edges. Let u £ V be an initial node, and 
c € Z<o be an initial credit. The decision problem of whether E(m) > c can be solved in 0(n ■ m) time and 0{n) 
space. 


5.2 The value problem for general graphs 

We now turn our attention to the value version of the minimum initial credit problem, where the task is to determine 
E(m) for every node u. The following claim establishes that if for all energies to reach some node v we have E„ {w) < 0, 
then E„(u) = —d{u, v), i.e., the energy to reach v from every node u is minus the distance from u to v. 

Claim 4. If for all w G V \ {w} we have Ev(,w) < 0, then for each u G V \ {t;} we have Ei,(?x) = —d{u, v). 

Proof. Let P ■. u vhs witness path to the distance, i.e., wt(P) = d{u, v) < oo (if d(u, v) = oo the statement 
is trivially true). Since every highest-energy node a; of P has E„(x) = 0, we have that x = v. Hence, P is a highest- 
energy prefix of itself, and for each prefix P' of P we have —wt(P) + wt(P') < 0 and thus E„(u) > — wt(P) = 
—d{u, v). By Remark]^ it is Ei,(m) < —d{u, v). The result follows. □ 

An 0{ri^ ■ m) time solution to the value problem. Claim|^together with Theorem|^lead to an 0{vf ■ m) method for 
solving the minimum initial credit value problem. First, we compute the X = {v GV : E(u) = 0} in 0(rif ■ m) 
time, by testing whether E(rt) > 0 for each node u. Afterwards, we contract the set A to a new node z, and by 
Claim|^for every remaining node u we have E(u) = E„(tt) = E^(m). Since u ^ X, the energy of u is 

strictly negative, and thus E 2 (u) < 0. Finally, by Claimj^ we have Ez(u) = —d{u, z). Hence it suffices to compute 
the distance of each node u to z, which can be obtained in O {n ■ m) time. 

In the remaining of this subsection we provide a refined solution of 0(k ■n-m) time, where fc = | A| +1 is the number 
of 0-energy nodes (plus one). Hence this solution is faster in graphs where k = o(n). This is achieved by algorithm 
ZeroEnergyNodes for computing the set X faster. 

Determining the 0-energy nodes. The first step for solving the minimum initial credit problem is determining the set 
X of all 0-energy nodes of G. To achieve this, we construct the graph G 2 = (V 2 , P 2 , wt 2 ) with a fresh node z 
follows: 

1. The node set is V 2 = 1^ U {z}, 

2. The edge set is P 2 = P U ({z} x V), 

3. The weight function wt 2 : E 2 Z is 


, . (0 \fu = z 

wt2(u, r’) = < / ^ , 

[ wt(tt, v) otherwise 

Remark 5. Since for every outgoing edge (z, x) of z we have wt 2 (z, x) = 0, if z is a highest-energy node in a path of 
G2, so is X. Hence every non-positive cycle in G2 has a highest-energy node other than z. 

Note that for every u G V, the energy E(u) is the same in G and G 2 . 

Algorithm ZeroEnergyNodes. Algorithm [^describes ZeroEnergyNodes for obtaining the set of all 0-energy nodes in 
G 2 . Informally, the algorithm performs a sequence of modifications on a graph g, initially identical to G 2 . In each 
step, the algorithm executes a Bellman-Ford computation on the current graph g with z as the source node, as long as 
a non-positive cycle G is discovered. For every such G, it determines a highest-energy node w of G, and modifies g 
by replacing every incoming edge {x, w) with an edge {x, z) of the same weight, and then removing w. See Figure]^ 
for an illustration. 

As 0-energy nodes are discovered, ZeroEnergyNodes performs a sequence of modifications to the graph g. We denote 
with g^ the graph g after the fc-th node has been added to X (and = G 2 ). We also use the superscript-A: in our 
graph notation to make it specific to g'^ (e.g. z) and E^(m) denote respectively the distance from u to z, and the 
energy to reach z from u in g'^). The following two lemmas establish the correctness of ZeroEnergyNodes. 
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Algorithm 3: ZeroEnergyNodes 

Input: A weighted graph G2 = (V2, -E2, wt2) 
Output: The set {w G V2 \ {z} : E(t)) = 0} 

1 Initialize sets V V2, E -B2 and map a/t t— wt2 

2 Let t? = (V, E, m) 

3 Initialize set X ■<— 0 

4 while True do 

5 Execute Bellman-Ford from source node z in 

6 if exists non-positive cycle C then 

7 Determine a highest-energy node w ^ z in C 

8 Assign X X U {m} 

9 foreach edge (x, m) G E do 
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if {x, z) ^ 'E then 

11 




Assign E ■<— E U {(*, z)} 

12 




Assign 'wt{x,z) t— wt2(x,w) 

13 



else 

14 




Assign z) t— min(wt2(a;, w), zyt(a;, z)) 

15 



end 

16 


end 


17 


Assign P •<— P \ {ui} 

18 

else 



19 


return X 

20 

end 



21 end 




Lemma 10. 

For every w G X we have E(t(;) = 0. 


Proof. The proof is by induction on the size of X. It is trivially true when |X| = 0. For the inductive step, let w be 
the k + 1-th node added in X. By line|^ m is a highest-energy node in a non-positive cycle C of We split into two 
cases. 


1. If z ^ C, then C is also a cycle of G, hence w is a highest-energy node in the infinite path V = of G, and 
E(u;) = 0. 

,z) = \Nt2{x,w') 


2. If z G C, let X be the node before z in C. By the modifications of lines 


11 


and 


14 


it is wGiyX, 


where w' is a node that has been added to X when the algorithm run on for some i < k.lt follows that w is a 
highest-energy node in a path P : z w' in G' 2 , and thus a highest-energy node in a suffix P' : w w' of P, 
where P' is a path in G. Hence Eiu/(u;) = 0. By the induction hypothesis, w' is a 0-energy node, i.e., E(r(;') = 0, 
thus by Claim|^we have E(u>) > Eiu/(r(;) = 0. 


The result follows. 


□ 


Lemma 11. For every w G V : E(u;) = 0 we have w G X. 

Proof. Consider any w G V : E(r(;) = 0. For some j G N, we say that “is aware of w” if either has a non-positive 
cycle C : w w, or w G X when |X| = i. Note that when ZeroEnergyNodes terminates there are no non-positive 
cycles in Hence, it suffices to argue that there exists a A: G N such that for each i> k, is aware of w. We first 
argue that there exists a k for which is aware of w. 

Let P be a witness for E{w) = 0, hence V traverses a non-positive cycle Ci in G, thus Gi exists in g^. Then there 
exists a smallest y G N such that some node w' of V is identified as a highest-energy node in a non-positive cycle 
C 2 (possibly Cl = C 2 ), and inserted to X. If w = w', we have that g^ is aware of w. Otherwise, since E(z(;) = 0 
and w' is a node in the witness V, we have Eu,/(zu) = 0. By the choice of w', the path V exists in gf therefore 
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Em'(w) = 0, and by Remarkwe have dP{w, w') < 0. It is straightforward that after the modifications 
and 14 we have that d^^^{w,z) < dP{w,w') < 0, and since wt^{z,w) = wt 2 (z,w) = 0, we have a 
non-positive cycle C : w w in through 0 . Hence either or is aware of w, thus there exists a /c £ N 
for which is aware of w. 


in lines 11 


Finally, observe that the distance z) does not increase in any for i > k until w is inserted to X, hence for 
each i > k, the graph is aware of w. The desired result follows. □ 


Determining the negative-energy nodes. Having computed the set X of all the 0-energy nodes of G, the second step 
for solving the minimum initial value credit problem is to determine the energy of every other node u G 1^ \ X. Recall 
the graph after the end of ZeroEnergyNodes. 


Lemma 12. For every u gV\X we have E(u) = {u, z). 


Proof. Consider any node u G V\X = \ {z}. By Cl aim| 4[ in the graph G we have E(it) = max^.E(„)=o Eu(u), 

and by the correctness of ZeroEnergyNodes from Lemma [Tofand Lemma 11 we have X = {v : E(ii) = 0}, thus 
E(m) = maxy^x E„(u). It is straightforward to verify that at the end of ZeroEnergyNodes, we have maxy^x Et,(ri) = 


EL^'(rt), i.e., the maximum energy to reach the set X in G is the energy to reach z in For all v G \ {z} 
it is eL'^'(ii) < 0, otherwise we would have E{v) = 0 and thus v G X and v ^ Then by Claim|^ eL'^'(m) = 
—(u, z). We conclude that E(m) = (u, z). □ 


Hence, to compute the energy E(u) of every node u G L \ X, it suffices to compute its distance to z in This 
is straightforward by reversing the edges of and performing a Bellman-Ford computation with z as the source 
node. Figure [^illustrates the algorithms on a small example. We obtain the following theorem. 



Fig. 3: Solving the value problem using operations on the graph Q. Initially we examine and a non-positive cycle 
is found (boldface edges) with highest-energy node x. Thus E(x) = 0, and we proceed with to discover E(u) = 0. 
In all cycles are positive, and the energy of each remaining node is minus its distance to z. 


Theorem 6. Let G = {V, E, wt) be a weighted graph of n nodes and m edges, and k = |{r; G V : E(r;) = 0}| -f 1. 
The minimum initial credit value problem for G can be solved in 0{k ■ n ■ m) time and 0{n) space. 


Proof. Lemma 10 Lemma[TT]and Lemma [T^ establish the correctness, so it remains to argue about the complexity. 
The while block of line [^ is executed at most once for each 0-energy node, hence at most k times. Inside the block, 
the execution of Bellman-Ford in line [^requires 0{n ■ m) time and 0(m) space. Since the Bellman-Ford algorithm 
uses backpointers to remember predecessors of nodes in distances, a highest-energy node w of a non-positive cycle 
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C in line|7]can be determined in 0{n). Finally, the/or loop of line|^will consider each edge {x, w) at most once, 
hence it requires 0{m) for all iterations of the while loop. Thus ZeroEnergyNodes uses 0{k ■ n ■ m) time and 0{n) 
space in total. The last execution of Bellman-Ford to determine the energy of negative-energy nodes does not affect 
the complexity. The result follows. □ 

Corollary 3. Let G = (V,E,wt) be a weighted graph of n nodes and m edges. The minimum initial credit value 
problem for G can be solved in 0{n^ ■ m) time and 0(n) space. 


5.3 The value problem for constant-treewidth graphs 

We now turn our attention to the minimum initial credit value problem for constant-treewidth graphs G = {V, E, wt). 
Note that in such graphs m = 0{n), thus Theoremgives an 0{vf) time solution as compared to the existing 
0(n‘* • log(n • W)) time solution. This section shows that we can do signihcantly better, namely reduce the time 
complexity to 0(n ■ logn). This is mainly achieved by algorithm ZeroEnergyNodesTW for computing the set X of 
0 -energy nodes fast in constant-treewidth graphs. 

Extended + and min operators. Recall the graph G 2 = (V 2 ,G 2 jWt 2 ) from the last section. Given Tree(G), a 
balanced and binary tree-decomposition Tree(G 2 ) of G 2 with width increased by 1 can be easily constructed by 
(i) inserting z to every bag of Tree(G), and (ii) adding a new root bag that contains only z. Let Z = Z x C x Z. For a 
map f : V 2 X V 2 —t Z, dehne the map p/ : C 2 x V 2 —> Z as 

u) < 0 or u = z 

\{f{u,v),v,f{u,v)) otherwise 

and for triplets of elements ai = (oi, 61 , ci), 02 = (^ 2 , ^ 2 , C 2 ) G Z, dehne the operations 

1 . min(ai,a 2 ) = ai with i = argmin^g^i 2 } aj 

2 . oi + 02 = («i + 02 : c), where c = max(ci, Oi + C 2 ) and 5 = 61 if c = Ci else h = & 2 - 

In words, if / is a weight function, then gf{u, v) selects the weight of the edge (m, v), and its highest-energy node (i.e., 
u if f{u, v) < 0 , and v otherwise, except when v = z), together with the weight to reach that highest energy node 
node from u. Recall that algorithm MinCycle from Section [^traverses a tree-decomposition bottom-up, and for each 
encountered bag B stores a map LDb such that LDb(u, v) is upper bounded by the weight of the shortest U-shaped 
simple path w -w u (or simple cycle, if u = v). Our algorithm ZeroEnergyNodesTW for determining all 0-energy nodes 
is similar, only that now LD^ stores triplets (a, b, c) where a is the weight of a U-shaped path P, & is a highest-energy 
node of P, and c the weight of a highest-energy prehx of P. For two triplets ai = (ai, 61 , ci), 02 = ( 02 , & 2 , C 2 ) € Z 
corresponding to U-shaped paths Pi and P 2 , min(Q;i,a 2 ) selects the path with the smallest weight, and ai + 02 
determines the weight, a highest-energy node, and the weight of a highest-energy prehx of the path Pi o P 2 (see 
Figure]^. 

Algorithm ZeroEnergyNodesTW. The algorithm ZeroEnergyNodesTW for computing the set of 0-energy nodes in 
constant-treewidth graphs follows the same principle as ZeroEnergyNodes for general graphs. It stores a map of edge 
weights wt : E 2 —i' ZU { 00 }, and initially wt{u,v) = wt 2 (u,u) for each {u,v) G E 2 . The algorithm performs a 
bottom-up pass, and computes in each bag the local distance map LD^ : B x B ^ I that captures U-shaped u v 
paths, together with their highest-energy nodes. When a non-positive cycle G is found in some bag B, the method 
KillCycle is called to modify the edges of a highest-energy node w of G and its incoming neighbors by updating the 
map wt. These updates generally affect the distances between the rest of the nodes in the graph, hence some local 
distance maps LDs need to be corrected. However, each such edge modihcation only affects the local distance map 
of bags that appear in a path from a bag B' to some ancestor B” of B'. Instead of restarting the computation as in 
ZeroEnergyNodes, the method Update is called to correct those local distance maps along the path B' -w B”. 

The following lemma establishes the correctness of ZeroEnergyNodesTW. Similarly as for Lemma[TO|and Lemma [TT| 
we denote with the graph obtained by considering the edges (u, v) for which wt{u, v) < 00 when |Ar| = k. 
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b pi 


Fig. 4: Illustration of the ai + operation, corresponding to concatenating paths Pi and P 2 . The path Pj denotes 
the z-th prefix of Pj. We have P = Pio P2, and the corresponding tripplet a = (a, b, c) denotes the weight a of P, its 
highest-energy node b, and the weight c of a highest-energy prefix. 


Algorithm 4: ZeroEnergyNodesTW 

Input: A weighted graph G 2 = (V 2 , E 2 ,wt 2 ) and a binary tree-decomposition Tree(G 2 ) 
Output: The set {« € V 2 \ {z} : E(u) = 0} 

// Initialization 

1 Assign X t— 0 

2 foreach u,v G V 2 do 

3 if («, v) e E 2 then 

4 I Assign w) t—wt 2 (M, u) 

5 else 

6 I Assign ■wt{u,v) t— 00 

7 end 

8 end 

// Computation 

9 Apply a post-order traversal on Tree(G), and examine each bag B with children Bi , B 2 

10 begin 

11 foreach u,v G B do 

12 I Assign LDB(M,t) t—min(LDs^(M,-!;), LDs 2 (M,ii),g„,(rt,u)) 

13 end 

14 if B is the root bag of a node x then 

15 foreach u,v G B do 

16 I Assign LD'b(m, t) t—min(LDs(it, t), LDB(ii, m) + LDB(a;, w)) 

17 end 

18 Assign LDb t— LDb 

19 if Bit e B with LDb (11, u) = (a, b, c) where a < 0 then 

20 Assign X t— X U {6} 

21 Execute KillCycle on b and B 

22 

23 

24 end 

25 return X 
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Methods: KillCycle 

Input: A 0 -energy node w and a bag B of Tree(G2) 

Output: Updates the local distance function LDs 

1 foreach edge {x,w) G E 2 do 

2 Assign 'wt{x, z) t— min(wt2(x, w), wtix, z)) 

3 Assign wt{x, w) t— 00 

4 Assign y t- argmax„g{2.,„} Lv(m) 

5 Let B' be the smallest-level ancestor of By examined by ZeroEnergyNodesTW so far 

6 Execute Update on By and its ancestor B' 

7 end 

8 return LDb 


Method 6: Update 

Input: A bag B' and an ancestor B" 

Output: The local distances LDb along the path B' ^ B" 


1 Traverse the path B' ^ B” bottom-up, and examine each bag B with children Bi, B2 

2 begin 

3 foreach u,v G B do 

4 I Assign LDb(u,u) min(LDsi(ti,-u), LDs2(M,f),p„,(M,ii)) 

5 end 

6 it B is the root bag of a node x then 

7 foreach u, u G B do 

8 I Assign LDb(m, u) t—min(LDB(ti, w), LDs(tt, *) +LDB(a:, w)) 


9 

to 

11 

12 

13 

14 

15 

16 end 


end 

Assign LDb -g- LDb 

if 3 m G B with LDs(m, m) = (a, b, c) where a < 0 then 
Assign X X U {b} 

Execute KillCycle on b and B 
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Lemma 13. For every v G V \ {z} we have v G X iff E(z;) = 0. 


Proof. We only need to argue that ZeroEnergyNodesTW correctly computes the non-positive cycles in every as 
then the correctness follows from the correctness Lemma [T0| and Lemma [TT]of ZeroEnergyNodes. Since by Remark[T] 
every cycle is a U-shaped path in some bag, it suffices to argue that whenever ZeroEnergyNodesTW examines a bag 
B (either directly, or through Update), every U-shaped simple cycle in B has been considered by the algorithm. This 
is true if no calls to KillCycle are made {if block in line 19 1 , as then ZeroEnergyNodesTW is the same as MinCycle, 
and hence it follows from Lemma 


Now consider that KillCycle is called and B' is the smallest-level bag examined by ZeroEnergyNodesTW so far. Let w 
be the 0-energy node, x an incoming neighbor of w, and y = argmax„g{2, Lv(u) (as in linej^of KillCycle). By the 
dehnition of U-shaped paths, the edge {x, w) appears only in paths that are U-shaped in bags along the path By B'. 
Hence, after setting wt{x,w) = oo (linej^of KillCycle), it suffices to update the local distance maps of these bags. 
Similarly, after setting wt{x, z) G- min(wt 2 (a;, w), wt{x, z)) (linej^of KillCycle), since B^ is the root of Tree(G' 2 ), it 
suffices to update the local distance maps in the bags along the path B^ B' . Either x = y, or, by the properties of 
tree-decompositions, B^ is an ancestor of By. Hence in either case B^ i?' is a subpath of By -w B', and both edge 
modihcations in linesandare handled correctly by calling Update on By and its ancestor B'. The result follows. 

□ 


Lemma 14. Algorithm ZeroEnergyNodesTW runs in 0{n ■ logn) time and 0{n) space. 


Proof. Let h = 0(log n) be the height of Tree(G 2 )- 


1 . 

2 . 

3. 


The method Update performs a constant number of operations to each bag in the path B' B" where B" is 
ancestor of B', hence each call to Update requires 0{h) time. 

The method KillCycle performs a constant number of operations locally and one call to Update for each incoming 
edge of w. Hence if w has kw incoming edges, KillCycle requires 0{h ■ kyf) time. Since KillCycle sets wt{x, w) = 
oo for all incoming edges of w, the node w will not appear in non-positive cycles thereafter. 

The algorithm ZeroEnergyNodesTW is similar to MinCycle which runs in 0{n) time and space (Lemma|^. The 
difference is in the additional if block in line 19 Since KillCycle is called when a non-positive cycle is detected, it 
will be called at most once for each node rt S V 2 \ {z} (from either ZeroEnergyNodesTW or Update). It follows 
that the total time of ZeroEnergyNodesTW is 


O 


J2ih-ku) 


0{n + h ■ \E 2 \) = 0{n ■ logn) 


where ku is the number of incoming edges of node u. Since KillCycle stores constant size of information in each 
bag of Tree(G 2 ), the 0(n) space bound follows. 


□ 


After the set X of 0-energy nodes has been computed, it remains to execute one instance of the single-source shortest 
path problem on the graph (similarly as for our solution on general graphs). It is known that single-source 
distances in tree-decompositions of constant treewidth can be computed in 0{n) time 0241191 . We thus obtain the 
following theorem. 


Theorem 7. Let G = (V, E, wt) be a weighted graph of n nodes with constant treewidth. The minimum initial credit 
value problem for G can be solved in 0{n ■ log n) time and 0{n) space. 


6 Experimental Results 

In the current section we report on preliminary experimental evaluation of our algorithms, and compare them to 
existing methods. Our algorithm for the minimum mean cycle problem provides improvement for constant-treewidth 
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graphs, and has thus been evaluated on low-treewidth graphs obtained from the control-flow graphs of programs. For 
the minimum initial credit problem, we have implemented our algorithm for arbitrary graphs, thus the benchmarks 
used in this case are general graphs (i.e., not constant-treewidth graphs). 


6.1 Minimum mean cycle 

We have implemented our approximation algorithm for the minimum mean cycle problem, and we let the algorithm run 
for as many iterations until a minimum mean cycle was discovered, instead of terminating after 0(log(n/e)) iterations 
required by Theorem]^ We have tested its performance in running time and space against six other minimum mean 
cycle algorithms from Tablej^in control-flow graphs of programs. The algorithms of Burns and Lawler solve the more 
general ratio cycle problem, and have been adapted to the mean cycle problem as in EtI . 



Madani 

Burns 

Lawler 071 

Dasdan-Gupta II26I 

Hartmann-Orlin 0^ 

Karp OSl 

Time 

O(n^) 

O(n^) 

0{n^ ■ log(n • W)) 

0{n^) 

O(n^) 

O(n^) 

Space 

0 (n) 

0{n) 

0{n) 

0{n^) 

O(n^) 

0(n^) 


Table 3: Asymptotic complexity of compared minimum mean cycle algorithms. 


Setup. The algorithms were executed on control-flow graphs of methods of programs from the DaCapo benchmark 
suit lO, obtained using the Soot framework ll45ll . For each benchmark we focused on graphs of at least 500 nodes. 
This supplied a set of medium sized graphs (between 500 and 1300 nodes), in which integer weights were assigned 
uniformly at random in the range { — 10^,..., 10^}. Memory usage was measured with m. 

Results. Figure 1^ shows the average time and space performance of the examined algorithms (bars that exceeded the 
maximum value in the y-axis have been truncated). Our algorithm has much smaller running time than each other 
algorithm, in almost all cases. In terms of space, our algorithm significantly outperforms all others, except for the 
algorithms of Lawler, Burns, and Madani. Both ours and these three algorithms have linear space complexity, but 
ours also suffers some constant factor overhead from the tree-decomposition (i.e., the same node generally appears in 
multiple bags). Note that the strong performance of these three algorithms in space is followed by poor performance 
in running time. 


6.2 Minimum initial credit 

We have implemented our algorithm for the minimum initial credit problem on general graphs and experimentally eval¬ 
uated its performance on a subset of benchmark weighted graphs from the DIMACS implementation challenges ifTl . 
Our algorithm was tested against the existing method of ifTOll . The direct implementation of the algorithm of ifTOl per¬ 
formed poorly, and for this we also implemented an optimized version (using techniques such as caching of intermedi¬ 
ate results and early loop termination). Note that we compare algorithms for general graphs, without the low-treewidth 
restriction. 

Setup. For each input graph we first computed its minimum mean value fi* using Karp’s algorithm, and then subtracted 
fj* from the weight of each edge to ensure that at least one non-positive cycle exists (thus the energies are finite). 

Results. Figure [^depicts the running time of the algorithm of ifTOl (with and without optimizations) vs our algorithm. 
A timeout was forced at 10^°/rs. Our algorithm is orders of magnitude faster, and scales better than the existing 
method. 
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Space (KB) Time (/js) Space (KB) Time (fts) 


2.4 

1.6 

0.8 

0.0 

1.2 

0.8 

0.4 

0.0 

2.4 

1.6 

0.8 

0.0 

1.2 

0.8 

0.4 

0.0 
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antir bloat chart eclipse fop hsqldb javac 


le3 



antir bloat chart eclipse fop hsqldb javac 


leS 



jflex jython luindex lusearch pmd xalan 



Fig. 5: Average performance of minimum mean cycle algorithms. 



Nodes 

Fig. 6: Comparison of running times for the minimum initial credit problem. 
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Madani 

Bums 

Lawler 

Dasdan-Gupta 

Hartmann-Orlin 

Karp 

Ours 

antlr 

55814 

61571 

165789 

284996 

21893 

7824 

18402 

bloat 

138416 

188356 

350302 

105145 

144171 

89949 

22391 

chart 

216962 

137112 

573767 

154062 

107229 

90717 

40890 

eclipse 

216859 

242323 

667869 

172792 

148523 

107864 

23486 

fop 

83080 

147384 

406371 

59176 

121742 

31557 

19306 

hsqldb 

131041 

153232 

208328 

86840 

228632 

40486 

19957 

javac 

58443 

110149 

122996 

179647 

14719 

34188 

20874 

jflex 

214297 

524822 

554093 

116820 

133323 

53329 

23860 

jython 

139106 

200922 

503766 

94052 

75569 

34864 

28760 

luindex 

199650 

217980 

1240411 

274319 

228856 

92379 

22142 

lusearch 

433211 

447280 

1180051 

263467 

333297 

101584 

55652 

pmd 

180551 

155118 

585315 

118578 

155682 

48326 

21978 

xalan 

120897 

156111 

394458 

81103 

96873 

47996 

14493 


Table 4; The time performance of Figure|^(in fj,s). 



Madani 

Bums 

Lawler 

Dasdan-Gupta 

Hartmann-Orlin 

Karp 

Ours 

antlr 

16805 

21018 

11144 

486435 

489176 

322384 

168648 

bloat 

29723 

24500 

19458 

1245272 

1249444 

826645 

306026 

chart 

27130 

30567 

18172 

2025448 

2029294 

1347048 

278586 

eclipse 

24215 

26488 

16293 

965063 

968595 

640720 

254393 

fop 

16845 

17975 

11052 

576174 

578646 

382338 

169738 

hsqldb 

16798 

19309 

11144 

486435 

489096 

322384 

168648 

javac 

14681 

17047 

9664 

372697 

375453 

247019 

144721 

jflex 

24561 

26946 

16322 

1244495 

1248036 

826743 

251549 

jython 

22518 

23337 

14899 

1059291 

1062570 

703581 

228207 

luindex 

39309 

40223 

25604 

3521607 

3526792 

2342833 

399076 

lusearch 

41488 

33350 

26991 

3387914 

3393343 

2253403 

422679 

pmd 

32204 

24481 

21021 

1391551 

1395786 

923975 

326137 

xalan 

16798 

17763 

11144 

486435 

489102 

322384 

168648 


Table 5: The space performance of Figure(in KB). 
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n 

Existing 

Existing Optimized 

Ours 

50 

9453565 

1680924 

48635 

58 

39744129 

3394193 

121774 

66 

55766874 

6201044 

267825 

74 

180080064 

12833610 

136239 

82 

267993314 

13563936 

116518 

90 

342779026 

25453589 

383292 

98 

74622910 

12648395 

501365 

106 

791441986 

60294150 

385799 

114 

1133055323 

80584700 

432290 

122 

1004898322 

67982455 

564838 

130 

2354354250 

165193753 

348112 

138 

881117317 

114743182 

636481 

146 

7050113907 

311146051 

501314 

162 

5179877563 

324877384 

1154447 

178 

Timeout 

589873640 

635155 

194 

3799301931 

391240954 

2672127 

218 

Timeout 

2596083382 

866213 

242 

Timeout 

2774469734 

1779512 

266 

Timeout 

2839496222 

7676638 

290 

Timeout 

6526762301 

1332403 

322 

Timeout 

5929433611 

1282258 


Table 6; The time performance of Figure in fj,s. 
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